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I.  INTRODUCTION 


The  coupling  of  electromagnetic  energy  from  one  region  to  another 
is  an  important  problem  in  many  areas  of  electromagnetic  engineering. 

Some  examples  are  leakage  from  microwave  ovens,  electromagnetic  penetration 
into  vehicles,  and  electromagnetic  pulse  interaction  with  shielded  electronic 
equipment. 

Our  approach  is  to  first  obtain  the  functional  equations  for  the 
general  aperture  to  cavity  to  aperture  coupling  system  using  the  equivalence 
principle  [1,  Sec. 3-5]  and  then  to  reduce  these  equations  to  a  matrix 
form  via  the  method  of  moments  [2].  The  various  matrices  are  interpreted 
in  terms  of  generalized  network  parameters,  such  as  voltages,  currents 
and  admittances  [2,3].  The  formulation  of  the  general  problem  is  similar 
to  the  formulation  of  Auckland  and  Harrington,  who  solved  for  coupling 
through  narrow  slots  in  thick  conducting  screens  [4]. 

Subsequently  the  problem  is  specialized  to  the  case  of  electrically 
small  circular  apertures  in  a  cylindrical  cavity  of  circular  cross  section 
aud  it  is  assumed  that  the  excitation  is  due  to  a  plane  wave  transverse  electromagnetic 
(TEM)  to  the  cylinder  axis.  The  aperture  admittances  are  obtained  using 
the  concept  of  polarizability  of  apertures  developed  by  Bethe  [5,6], 
and  the  radiation  terms  introduced  by  Harrington  [7].  In  general  an  electrically 
small  aperture  can  be  described  in  terms  of  a  moment  solution,  where 
at  least  three  expansion  functions  must  be  used  to  express  the 
equivalent  magnetic  current.  However,  for  simplicity,  we  assume  that 
symmetry  is  such  that  one  expansion  function  will  suffice. 

The  cavity  region  can  be  viewed  as  a  short  circuited  cylindrical 
waveguide  [8].  Thus  waveguide  theory  can  be  applied  and  the  field  in  this 
region  expanded  in  terms  of  circular  waveguide  modes.  Furthermore, 


we  assume  that  all  waveguide  modes  except  the  dominant  TE^  mode  are  in  the 
cutoff  condition,  thereby  obtaining  relatively  simple  formulas  and  a  simple 
equivalent  circuit  for  the  coupling  problem. 

Special  attention  has  been  paid  to  the  case  of  TE  (transverse  electric 
to  the  cylinder  axis)  oblique  incidence  upon  the  structure  and  to  the 
effects  of  cavity  losses,  both  lossy  material  filling  the  cavity  and 
finite  conductivity  of  the  cylindrical  conductor,  on  the  transmission 
coefficient. 

Finally,  our  discussion  is  extended  and  we  assume  that  two  waveguide 
modes,  namely, the  TE^  and  the  TM^,  propagate.  TMQ1,  TE^.and  TEQ1  modes 
can  propagate  but  are  not  excited. 


II.  FORMULATION  OF  THE  PROBLEM 

The  problem  to  be  considered  is  that  of  coupling  between  two  half-space 

regions  through  a  cylindrical  cavity.  Pig.  1  illustrates  the  cross  section 

of  the  geometry  of  the  problem.  The  left-hand  half  space  (z  <  0)  is 

called  region  a,  the  circular  cavity  region  (0  <  p  <  a,  0  <  z  <  d) 

is  called  region  b,  and  the  right-hand  half  space  (z  >  d)  is  called  region 

c.  A  circular  cavity  is  a  cylindrical  cavity  of  circular  cross  section. 

The  boundary  common  to  regions  a  and  b  is  called  the  aperture  A^.  The 

boundary  common  to  regions  b  and  c  is  called  the  aperture  A Regions 

a,b,  and  c  are  each  filled  with  homogeneous  media  of  constitutive  parameters 

(y  ,  e  ) ,  (y, ,  £,),and  (y  ,  e  ),  respectively.  We  are  not  considering 
a  a  b  d  c  c 

dissipation  and  therefore  each  y  and  each  c  Is  real.  The  excitation  is  due  to 
known  sources  J_^  and  M*",  with  exp(jait)  time  dependence,  in  region  a. 

The  equivalence  principle  [1,  Sec.  3-5]  is  used  to  divide  the  original 
problem  into  three  equivalent  problems,  as  shown  in  Figs.  2-4.  In  region 
a,  the  field  is  produced  by  the  original  sources  J_i,  M*,  plus  the  equivalent 
magnetic  current  -M ^ ,  where 

M,  -  n  x  E  (1) 

—1  —a  — 

over  the  aperture  region  A^ ,  all  radiating  in  the  presence  of  a  complete 
conductor  (aperture  A^  shorted).  In  (1),  n^  is  a  unit  vector  normal  to 
A^.  and  E  is  the  electric  field  in  the  aperture  A^  in  the  original  problem. 

In  region  b,  the  field  is  produced  by  the  equivalent  currents  Mj,  given  by 
(1)  over  A-p  and  -Mj,  where 

M  *  E  *  n  (2) 

— 2  —  — c 

over  the  aperture  region  A2,  both  radiating  in  the  presence  of  a  conductor 

completely  enclosing  the  cylindrical  cavity  region  b  (both  apertures 

shorted).  In  (2),  the  unit  vector  is  normal  to  Ay,  and  E  is  the  electric 

field  in  the  aperture  A2  in  the  original  problem.  Finally,  in  region  c  we 
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z  »  0  z  *  d 


?ig.  1.  Geometry  of  the  coupling  problem  between  two  half-space 
regions  through  a  circular  cavity. 
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Fig.  2.  Equivalence  for  region  a. 
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have  the  equivalent  magnetic  current  given  by  (2)  over  the  aperture 
region  A^,  radiating  in  the  presence  of  a  complete  conducting  plane 
(aperture  shorted) . 

The  use  of  -M^  in  region  a  and  in  region  b  ensures  continuity  of 
the  tangential  components  of  the  electric  field  across  the  aperture  Aj . 

The  use  of  -M?  in  region  b  and  M?  in  region  c  ensures  continuity  of  the 
tangential  components  of  the  electric  field  across  the  aperture 
Continuity  of  the  tangent ial  components  of  H  across  each  aperture  leads  to 
the  operator  equations  for  the  problem.  The  procedure  is  described  in  detail 
in  [3].  The  result  is 

-H®(Mi)  -  H^Mj)  +Ht%)  =  -HjC  over  A^ 


-  h^(m2)  -  h^(m2)  =  o 


(3) 


over  A2 


where  HP(M  )  denotes  the  tangential  component  of  H  due  to  M  radiating  in 
— t  — q  —  — q 


,sc 


region  p  with  all  apertures  shorted,  and  Ht  denotes  the  tangential  component 
of  H  due  to  the  impressed  sources  J* ,  M*  in  region  a  with  aperture  A^  shorted. 
Equation  (3)  is  first  solved  for  the  equivalent  magnetic  currents  and  M2> 
and  then  the  fields  in  each  region  can  be  computed  from  these  equivalent 
currents. 

To  obtain  a  solution  to  (3),  it  is  convenient  to  use  the  method  of 

moments  [2].  The  procedure  is  described  in  detail  in  [3].  The  result  can 

be  summarized  as  follows:  Sets  of  expansion  functions  {M^}  in  A^  and 

{M_  }  in  An  are  assumed,  and  the  equivalent  currents  are  expressed  as 
— 2n  2 


M.  =  \  V.  ML 
—1  **,  In— In 

n=l 


s2 

—2  •  l  'A 

n=l 


(M 


8 


and  substituted  into  (3).  Symmetric  products 
<  A,  B  =  |  |  A  •  B 


ds 


(5) 


<  C,  D  >2 


C  •  D  ds 


are  defined  for  each  aperture.  Sets  of  testing  functions  {W^}  in  A^  and 

{W„  }  in  A0  are  defined,  and  equations  (3)  are  tested  with  each  W  , 

— /n  z  — qtn 

m  =  1,2,. ...N  .  The  result  is 

q 


where 


<Yu>  Yi  +  'Yn>  ?1  +  'Yn'  Y2  ■  fl 


!Y21J  Y1  +  '41  f2  +  'Y22>  ?2  *  3 


[YP  ]  =  [-  <w  ,  HP(M  )>  ] 
qq  -t  -qn  qJN  X  N 


(6) 


[Y„ J  =  t  <w  »  HP(M  )>  1 
qr  -qm  —  t'-rn7  q  x  fj^,  q  f  r 


(7) 


(8) 


5,  '  <Vn„  x  1 

q 


(9) 


The  matrices  [YPr]  are  called  the  generalized  admittances,  the  vector  I*  is 
called  the  generalized  source  current,  and  the  vectors  are  called  the 
generalized  voltages.  A  solution  of  the  problem  is  obtained  by  solving  the 
matrix  equations  (6)  for  and  which  determine  the  equivalent 
magnetic  currents  and  by  (A).  If  a  Galerkin  solution  is  used,  that 

is,  if  {W^}  *  {M^n)  an<*  2n^  =  2n^  it:  t*len  follows  from  (7)  and  (8)  that 


9 


t 


(10) 


[<,1  •  [-  <v  a?*",,)  >  * „o 


[Y**  ]  =  [  <M  ,  H^(M  )  >  Jx,  x  H  , 

<lr  — qm  —  t  —  rn  <T  N  x  N  ,  q  f  r 


■  i  -  <«lm.  «r  > !  jKj  x ! 


(id 


It  is  important  to  note  that  computation  of  [Y^]  involves  only 

|j  I  i 

region  a  ,  computation  of  (Y^J,  [Y^J,  ^Y21^,and  [Y^]  invo*ves  only  region  b, 
and  computation  of  [Y22]  involves  only  region  c.  Hence,  we  have  divided 
the  problem  into  three  parts,  each  of  which  may  be  formulated  independently, 
and  therefore  we  can  use  previous  computations  of  problems  having  one 
of  these  parts  in  common.  For  example, the  aperture  admittances  of  electrically 
small  holes  for  radiation  into  half-space  will  be  used  in  the  computation  of 
[Yjj]  and  [Y^l  * 
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reduces  to  the  single  element 


Ait 


11  2>yaaml  3n  X2 


(15) 


a  a 


where  a  ^  is  the  y-directed  magnetic  polarizability  of  a  circular 
aperture  related  to  its  radius  by  [9] 


a  .  =  4ri 

ml  3  1 


(16) 


and  p  ,  n  ,  and  X  are  the  permeability,  intrinsic  impedance,  and  wavelength 
el  3l  3 

of  medium  a,  respectively. 

In  region  b,  the  cavity  part  of  the  problem, excites  the  TE  and 

TM  circular  waveguide  modes  having  a  modal  magnetic  vector  with  nonzero 

y  component  at  the  cylinder  axis.  The  modes  excited  by  are  the  TE^p 

and  TM.  for  all  P.  These  modes  will  excite  only  the  y-directed 
ip 

magnetic  dipole  mode  of  the  aperture  k^.  Hence  we  let 


~2  =  V21«21 


(17) 


where  is  the  quasi-static  current  which  produces  the  effect  of 
a  unit  magnetic  dipole  K£  =  1  in  the  y  direction,  and  is  a 
complex  coefficient  to  be  determined. 

The  transverse  (to  z)  components  of  the  electric  and  magnetic  fields 
due  to  are  called  (M^)  and  H^M^) ,  respectively,  and  are  expanded 
in  terms  of  waveguide  modes  as  [1,  Sec.  8-1]  and  [8] 


it«rl>  -  l  VrU(z)  ^ 

V4l>  ■  {  'nS*  i 1 


(18) 


(19) 
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Here  e^  and  are  the  normalized  modal  electric  and  magnetic  field  vectors 
and  Vf^(z)  and  are  c^e  moc*al  voltages  and  currents  due  to  only, 

calculated  for  the  equivalent  problem.  The  summations  in  (18)  and  (19) 
are  assumed  to  be  over  all  modes,  both  TE  and  TM  if  necessary.  The  total 
transverse  (to  z)  electric  and  magnetic  fields  are  then  given,  respectively , by 


St  •  ''uW  -  V!A%1> 


(20) 


and 


«t  •  »u«t®ll>  -  V2lV^21> 


(21) 


In  the  one-term  Galerkin  solution,  the  cavity  matrix  of  admittances 
becomes 


yb 

11  12 


yb  yk 

21  *22 


I-ivlli(0)  Y.  cot  k±d 


I-JViii(O)  V21i(d)  Y1  CSC  ktd 


?"jVlli(0)  V21i(d)  Yi  CSC  kid 


Hv2ii(d)  Yi  cot  kid 

(22) 


where  Yi  is  the  ith  waveguide  mode  characteristic  admittance,  k±  is  the 
ith  waveguide  mode  number,  and 


vm(°)  = 


apert 

A1 


h .  ds 
— i 


(23) 
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i«>  ■  \ 


—21  *  — i  ds 


apert 

A„ 


(24) 


If  all  waveguide  modes  except  the  dominant  TE^  (i  =  1)  mode  are  in  the  cutoff 
condition,  then  all  Y^  and  ^  are  imaginary  except  and  k^.  Note  that 
there  are  two  TE.^  modes,  but  only  one  of  them  is  excited.  We  let 


\ 


14 


Yi-  ^Bi 


ki  -  -jai 


i  t  1 


(25) 


and  the  elements  of  the  admittance  matrix  in  (22)  become 


Yb  Yb 

11  12 


Yb  Yb 

Y21  22  J 


-Jn^lYl  cot  k^d  +  jB11 
jnlln21Yl  CSC  kld  +  3B21 


JnHn21Yl  csc  kld  +  3B12 

”^n21Yl  COt  kld  +  ^B22 

(26) 


where 


nll  =  Vlll(0)’  n21  V211(d) 


(27) 


'll 


t  vi i 4 (°)®4  coth  ad, 
i#l 


22 


i  v: 

i*l 


21i(d)Bi  coth  aid 


(28) 


B12  B21 


Jl 


vni(0)  V„„(d)  csch  c^d 


221 


(29) 


Whenever  >>  1  (i  4  l).csch  a^d  -*•  0.  Hence,  B^  and  B,^  may  be 
neglected  in  Yb^  and  Yb^,  respectively.  This  fact  is  evident  since  the  fields 
due  to  those  nonpropagating  modes  exist  primarily  in  the  vicinity  of  the 
apertures,  and  thus  their  contribution  to  the  coupling  between  the  two 
apertures  is  negligible.  However,  these  fields  do  contribute  to  the  input 
admittances  Yb^  and  Yb2  of  the  two  cavity  backed  apertures.  As  the  sizes 
of  the  apertures  decrease, the  fields  become  more  localized,  more  modes 
are  required  for  an  adequate  respresentation ,  and  the  proportionate 
contribution  of  each  particular  mode  (such  as  the  1=1  mode)  becomes 


\ 


'  ^ 
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smaller.  Therefore,  provided  the  cavity  sides  are  not  close  to  other 
boundary  surfaces,  the  susceptances  B^  and  B22  are  essentially  the 
same  as  the  susceptance  of  an  aperture  in  an  infinite  plane  conducting 
screen  [10].  Hence  we  let. 


3  ’ 

8wy,  R, 

b  l 


8u)ybR2 


=  B2i  =  ° 


where  y,  is  the  permeability  of  medium  b.  The  matrix  of  admittances 

D 

Eq.  (26)  becomes 


2 

~^nllYlCOt  kld  +  ^Bll 
^nlln21Yl  CSC  kld 


•^nlln21Yl  CSC  kld 

2 

-jn2iYi  cot  kjd  +  jB22 


(31) 


This  is  the  two-port  admittance  matrix  for  a  transmission  line  of 


length  d,  having  a  characteristic  admittance  Y^,  and  a  propagation 
constant  k^,  with  a  shunt  susceptance  connected  to  its  input  terminal 
through  an  ideal  transformer  of  turns  ratio  n^  and  a  shunt  suspeptance 
B22  connected  to  its  output  terminal  through  an  ideal  transformer  of 
turns  ratio  This  two-port  matrix  is  the  part  of  the  equivalent  circuit 

labeled  region  b  in  Fig.  6. 

The  characterization  of  region  c  is  similar  to  that  of  region  a. 

Hence ,  one  has 


3  +  4tt 

j8wy  r3  3t)  A2 
c  2  c  c 


(32) 


where  yc»  nc,  and  Ac  are  the  permeability,  intrinsic  impedance,  and  wave¬ 
length  of  medium  c,  respectively. 

Thanks  to  (31),  the  equivalent  circuit  for  the  one  term  moment  equation  (6) 
is  the  circuit  given  in  Fig.  6.  This  circuit  may  be  reduced  to  the  equivalent 
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circuit  of  Fig.  7  because 


Mn  •  hx  ds 


-11  *  Vs 


"2  "1 

so  that  for  practical  considerations  we  may  take 


n21  nll 

A  parameter  of  interest  is  the  transfer  admittance  [4] 

x  yB  yB  (Ya  +  yB  )  (YB  +  Yc  ) 

I  12  21  U11  11M  22  22; 

12  V  b  Kii) 

21 

which  allows  one  to  calculate  the  strength  of  M2  =  V21— 21’  8iven  t^ie 
excitation  I1  of  (14).  Setting  n21  =  n^  in  (31)  and  then  substituting 
(31)  into  (33),  we  obtain 

Y12  =  [Y11  +  Y22  +  j(Bll  +  B22)1cos  kld 


(34) 


2..  ,  <Yn  +  %l><Y22  +  jB22>  .  .  , 

+  ^  nilYl  +  - 2 - Jsin 

nnYi 

Equation  (34)  can  also  be  obtained  from  the  equivalent  circuit  in  Fig.  7. 
Denoting  the  aperture  admittances  by 

Y11  =  ^  +  j®3  and  Y22  =  gC  +  J®C  (35> 

we  can  write  the  real  and  imaginary  parts  of  Y12  as  given  by  (34)  as 

Ga(BC  +  B  )  +  GC(B3  +  B  ) 

Re(Y12)  =  (Gd  +  G  )cos  - - - —  sin  kd  (36) 

nllYl 


Im(Y12)  =  (Ba  +  BC  +  Bn  +  B22)cos  k  d 


,  GaGC  -  (B3  +  B  )(BC  +  B„) 

+ - rr1 - — V 

nllYl 


(37) 
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To  complete  our  formulation  we  need  to  specify  the  parameters  re¬ 
lated  to  the  i  =  1  (TE^)  dominant  mode  of  the  circular  waveguide. 

The  normalized  modal  magnetic  and  electric  field  vectors  of  the  TE^  mode 
in  a  cylindrical  coordinate  system  are  [1,  Sec.  81] 


„  o o  ~r  r  1.841  T '  ,  1 . 841  ,  ,  ^  .  1  T  ,1.841  ,  ^  , 

h.i  =  -0-887[  — - —  JL( — - —  P)  sin  4>  — —  p)  cos  <p 

e,  =  h,  x  u 
— 1  —1  —  z 


(38) 

(39) 


where  J^(x)  is  the  Bessel  function  of  the  first  kind  of  order  1  and  J^(x) 
denotes  the  derivative  of  J^x)  with  respect  to  x.  The  modal  wave  number  is 


r  2  ,1.841,2.1/2 

kl  =  ^  Ubeb  -  ] 


(40) 


and  its  modal  characteristic  admittance  is 


Y  =  — - 
1  wy. 


(41) 


Another  parameter  of  interest  is  the  turns  ratio  n^.  Substitution  of  (38) 
into  (23)  yields 

0.816 


vm(0)  *  . 


(42) 


and  from  (27)  one  readily  obtains 


0.816 


lll 


(43) 
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IV.  POWER  TRANSMITTED 


The  power  transmitted  through  the  cavity  to  region  c  is  equal  to 

c 

the  power  dissipated  in  Y 22  of  the  equivalent  circuit,  that  is 


Ptrans  ’  lV2ll2  G' 

In  terms  of  the  transfer  admittance  (33),  this  is 


(44) 


trans 


r1  2 

lf-1  GC 

12 


(45) 


where  I*  is  given  by  (14),  Y^  by  (36)  and  (37),  and  Gc  =  Re(Y22).  Here, 


Y22  is  given  by  (32), 


The  power  incident  upon  the  aperture  when  the  incidence  is 


normal  is 


p”  =  n  Jh  | 

me  a 1  o' 


(46) 


where  Hq  is  the  amplitude  of  the  incident  field,  and  A  =  ttR^  is 
the  area  of  the  aperture  A^.  The  superscript  n  means  that  the  power 
defined  here  is  for  normal  incidence.  One  can  define  the  transmission 
coefficient  T  of  the  cavity  to  be  the  power  transmitted  through  the 
cavity  to  region  c  normalized  witli  respect  to  the  incident  power  (46) , 
that  is 


T  = 


P 

trans 


(47) 
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Now,  substitution  from  (14),  (45)  and  (46)  into  (47)  results  in 


na  I Y12 1  A  3Vc  Ac  lYl2|2  A 


The  transmission  coefficient  (48)  depends  on  the  cavity  depth  d 
only  through  the  parameter  Y^.  It  will  attain  its  maximum  when  |Yl2| 
is  minimum.  For  small  apertures  we  see  that  Ba>>G3  and  BC>>GC.  Hence,  the 
coefficients  of  the  trigonometric  terms  in  (37)  are  much  larger  than  those 
in  (36),  and  we  can  minimize  [yi2 |  by  setting  lm(Yl2)  *  0.  Doing  this 
and  retaining  only  dominant  terms,  we  obtain 

(B3  +  BC  +  B  +  B  ) 

tan  k,d  = - — - —  n  y  (49 

1  a  c  ill 

(B  +Bn)(B  +  B22)  111 

Since  the  right-hand  side  of  Eq.  (49)  is  a  small  negative  number,  a 
first  order  approximation  to  resonance  is 


(B3  +  BC  +  B„  +  B„) 

k.  d  -  mu  +  - — - —  n  Y 

1  (Ba+Bll)(BC+B22)  111 


„  .  ^  Bn  *  Btt>  jiL, 

”S  (Ba  +  Bu)(.'  +  B22)  *  2 


Here  the  subscript  "res"  denotes  "at  resonance",  A  is  the  wavelength 
of  the  dominant  mode  in  the  cavity  region,  and  m  =  1,  2,  3 . 
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Next  we  wish  to  obtain  the  transmission  coefficient  at  resonance 


Using  the  approximations 


c°s  k1dreg 


-  (-1) 


(51) 


sin 


k1  d 
1  res 


(-l)m 


(B  +  B  +  B. .  +  0 

- 11 - n“Y 

(B3  +  BU)(BC  +  B22) 


in  (36),  we  obtain 


12  res 


Re (Y  _)  ,  =  (-I)1"  ^  (G3  +  GC)'  -  (Ba  +  BC  +  Bn  +  B22)[ - 5 -  + 


11  22  Ba  +  Bu  BC  +  B 


If  we  introduce  the  parameters 

-  _  ^ 
^b  M 

a 


Cb  "  r 


(52) 


(53) 


yc  =  IT 


rc  =  F 


R  ^ 

(R  } 
R2 


(54) 


we  obtain 


,Vi/2 

uc  =  ir) 

c 


(55) 


X  = 


(ye) 

c  c 


1/2  a 


(56) 
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Now  substituting  for  Ga,  B3,  B^,  B ^ *  G°  and  8°  as  obtained  from  (15), 

(30)  and  (32),  using  the  parameters  defined  in  (53)  and  (54),  and  applying  the 
relations  (55)  and  (56),  we  cast  (52)  into  the  form 


m  4  it 


2 

3n  x 
a  a 


1  +  €  (y  £  ) 

c  c  c 


1/2 


-(X(hb  +  Mc)  +  hc(l  +  hb) ] [ 


e  (5  c-  )1/2 

c  c  c 


Uc(l  +  yb)  x(yb  +  uc) 


(57) 


At  resonance,  (48)  can  be  written  as 


16tt  £  (jl  e  )1^2 
c  c  c 

rreS  =  3nJ  A2  !Re(Y12)res|2  A 


(58) 


and  (50)  as 


X(U,  +  U  )  +  U  (1  +  y  )  R,  X 
^  ,  lm  .  3.556  -JL_]  -j- 


(l  +  yb)(ub  +  yc) 


a  A, 


(59) 


If  medium  a, medium  b,  and  medium  c  are  the  same  (i.e.,  y,  =  y  =  £ 

b  c  c 


1) 


and  if  Rj  =  R2  (i.e,,  X  =  1)  >  (57)  reduces  to 


Re 


nri-1  87r 


3n  X2 
a  a 


(60) 


The  resonance  occurs  at 


res 


[m  -3.556 


Ri  A 

v2T 


(61) 


and  the  transmission  coefficient  at  resonance  (58)  becomes 


3X 

r  =  —a 

res  4ttA 


(62) 
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The  transmission  cross  section  of  the  system  is  defined  as  the 
area  for  which  the  incident  wave  contains  the  power  transmitted  by  the 
system.  It  follows  that  the  transmission  cross  section  is  equal  to  TA. 
From  (48) ,  one  obtains 


"a^nl2  !'A  xclyi2l2 


If  media  a,  b,  and  c  are  the  same,  one  obtains  from  (62)  that  the  trans- 

2 

mission  cross  section  at  resonance  becomes  (TA)  =  3 A  / 4tt  independent 

res  a 

of  the  common  radius  R^  =  of  the  apertures.  This  transmission  cross 
section  and  the  transmission  coefficient  (62)  are  the  same  as  those  obtained 
in  [7]  for  a  small  slot  in  a  zero-thickness  plane  conducting  screen, 
resonated  by  a  capacitor  placed  across  its  midpoint,  and  excited  by  an 
incident  electric  field  perpendicular  to  its  axis. 


OBLIQUE  INCIDENCE  UPON  THE  STRUCTURE  (TE  CASE), AND  ONK  PROPAGATING 


MODE  IN  THE  CAVITY 


Consider  a  plane  wave  incident  upon  the  structure  at  some  angle 
in  the  y-z  plane  measured  from  the  negative  z  axis,  as  shown  in 
Fig.  8.  This  is  a  TE (transverse  electric  to  the  cylinder  axis)  excitation. 
Following  an  approach  similar  to  that  of  section  ITT,  we  reduce  the  excitation 
vector,  Eq.  (11),  to  the  scalar 

I  =  -H^C(0)  =  -2H  cos  0 

y  o  inc 


where  Hq  is  the  amplitude  of  the  incident  magnetic  field. 

Except  for  the  fact  that  I1  is  nowglvenbv  (64)  instead  of  (14),  the 
equivalent  circuit  and  its  parameters  remain  unchanged.  The  power  transmitted 
through  the  cavity  to  region  c  is  given  by  (44) ,  and  in  terms  of  the 
transfer  admittance  bv  (45).  Substitution  of  I1  from  (64)  into  (45)  results  in 


P 

trans 


G 


c 


where  Y^  is  given  by  (36)  and  (37)  and  G°  =  lUKY^)  with  Y^  given  by 

The  transmission  coefficient  is  defined  by  Eq.  (47).  Substituting 
from  (65)  and  (46)  into  (47) ,  one  obtains 


(65) 


(32). 


T  » 


4G 


2 

—  cos  0. 

2  inc 


16tt 


\  'yi21  a 


3naVclYi2l  A 


cos^  9 


inc  (66) 


Note  that  here  ,  the  actual  power  incident  upon  the  aperture  A^,  is 


related  to  the  power  p^nc  given  by  (46)  by 


VI.  THE  EFFECT  OF  CAVITY  LOSSES 


conductivity.  Equating  the  attenuation  constants  of  (69)  and  (70), 
we  obtain  that  a  lossy  conductor  case  behaves  like  a  lossy 

I  II 

dielectric  case  (e^  =  -  je^  and  quality  factor  =  r^/e^')  if  the 

quality  factor  due  to  conductor  losses,  defined  as. 


kb  nb  a 
2R 


t(^)2  + 

V 


-] 


(1.841)*  -  1 


is  equal  to  the  quality  factor  due  to  dielectric  losses  Q, . 


VII.  TOO  PROPAGATING  MODES  IN  THE  CAVITY  (TEM  AND  TE  EXCITATIONS) 


The  problems  to  be  considered  are  the  problems  of  Fig.  5  (TEM 
excitation)  and  of  Fig.  8  (TE  excitation)  treated  in  sections  III  and 
V,  respectively.  Here  we  extend  our  discussion  and  assume  that  two  excited 
waveguide  modes,  namely,  the  TE.^  and  the  TM^,  propagate.  The 
differences  appear  in  the  admittance  representation  of  the  cavity  region. 

The  parameters  of  the  equivalent  circuits  for  region  a  and  region  c 
remain  unchanged,  and  may  be  extracted  from  sections  III  and  V  for  the 
TEM  and  TE  excitations,  respectively. 

If  all  excited  waveguide  modes  except  the  TE^  (i  =  1)  and  TM^  (i  =  2) 
modes  are  in  ti.e  cutoff  condition,  then  all  Y^  and  k^  are  imaginary  except  , 
k,  ,  Y„,  and  k_.  We  let 

12  2  v 


Yi  "  ^Bi 


ki  -  -j«t 


i  4  1,2 


and  the  elements  of  the  admittance  matrix  in  (22)  become 


(72) 


Yb  Yb 
11  12 


"iVn1(°)Y1  cot  k:d 

-JV^12(0)y2  cot  k2d+;*Bil 


-jVlll(0)V211(d)YlCSC  kld 
-jV112(0)V212(d)Y2  C8C  k2d+jB12 


'jVlll(0)V211(d)Yl  CSC 

k^d 

-jV2U(d)YlCot  kjd 

-jVn2(0)V212(d)Y2  esc 

k^+jS’i 

-jV2212(d)Y2  cot  k2d+jB^2 

(73) 
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where 


=  l  V?. . (0) B .  coth  ad; 

tfl,2  111  1 


B22 


=  I  V*  (d)B  coth  ad 

i>l,2 


(74) 


R'l2  =  B21  2  Vlli(0)  V21i(d)Bi  CSCh  V 

a  .  b  ,  b 

Whenever  a  d»l  (i  *  1,2),B|2  and  may  be  neglected  m  Y  n  andYn> 

respectively.  Using  an  argument  similar  to  that  of  section  III,  we 
conclude  again  that  B^  and  B^  are  essentially  the  same  as  the  sus- 
ceptance  of  an  aperture  in  an  infinite  plane  conducting  screen.  Thus 
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11 


3  o  t  _  _ 

3  ’  22 


8wVbR1 


3 

8u)MbR2 


(76) 


Now,  the  matrix  of  admittances  (73)  becomes 

2  2 
4  ^  V  1.  A  2  _  V  1,  A  , 


r  b 

„b  ‘ 

Y 

11 

Y 

12 

„b 

„b 

Y 

21 

Y 

22 

Jnlln21Yl  «c  k1d+jn12n22Y2  esc  k2d 


-jn^cot  kbd -jn22Y2  cot  k.,d +jB’2 

(77) 


where 


nll  =  Vlll(0) 


n12  V112(0) 


n21  =  -V2U(d) 


( 7  f  ) 


"22  =  -V212(d> 


This  is  the  two- port  admittance  matrix  for  two  transmission  lines  of  length  d, 
having  characteristic  admittances  Y^  and  Y, , respect ively. and  propagation 
constants  k^  and  k2 .respectively ,  with  shunt  susceptances  Bjj  and  B22 


connected  through  ideal  transformers,  at  their  input  and  output  ports. 

The  equivalent  circuit  for  this  one  term  moment  solution  is  given 
in  Fig.  9,  and  it  may  be  reduced  to  the  equivalent  circuit  of  Fig.  10. 


This  last  step  is  due  to  the  fact  that  for  the  propagating  modes  we 
have 


-11 


ds 


M 2i*  h i  ds 


and 


Mn  .  h2  ds 


—21"  —2  ds‘ 


Therefore,  for  practical  considerations,  we  may  take  n .  =  n  .  and  n  =  n 

21  11  22  12 

The  transfer  admittance,  Y^2>  defined  by  Eq.  (33)  is  obtained  from 
Fig.  10.  That  is. 


=  yYii  +  Y22  +  J(Bli+  ®>2^nllYlC0S  •'id  sin  k2d  +  ni2Y2  cos  k2d  sin  kid  ] 


+  j[(nnV2  +  Cn12Y2) 2  +  (Y11  +  jBll)(Y22  +  JB22)]sin  kld  sin  k2d 


.2  2  1  2  2  -1 
+  2Jnnni2YlY2  [1  -  cos  k^d  cos  k,d)>  [n^Y^  sin  k2d  +  n^Y,,  sin  kid] 

(79) 


Once  the  transfer  admittance  is  obtained  from  (79)  the  transmission 
coefficient  may  be  computed  from  (48)  and  (66)  for  normal  and  oblique 
incidence,  respectively. 
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The  parameters  related  to  the  TE^  (1  =  1)  and  TM^  (1=2)  modes  should 

he  specified.  The  normalized  magnetic  and  electric  field  vectors  of  the 

TE^  mode  in  a  cylindrical  coordinate  system  ,  as  given  by  (38)  and  (39), 
are 


hj  =  -0.887[^^  .l|(j^iP)  sin  (D  up  +  ^  cos  4>  ^  (80) 


e,  =  h,  *  u 
—1  ~1  — z 


Its  modal  wave  number,  as  given  by  (40),  is 


r  2  ,1.841x2,1/2 

kl  =  Vb  "  (_7_)  1 


and  its  modal  characteristic  admittance,  as  given  by  (41) ,  is 


1 


The  normalized  magnetic  and  electric  field  vectors  of  the  TM^  mode 


in  a  cylindrical  coodinate  system  are  [1,  Sec,  8-1] 


„  r,„,l  ,  ,3.832  x  .  .  .  3.832  ,.,3.832  , 

h  =  -0.517  [—  J  (- — - —  p)  sin  <p  u  +  — - —  J'  ( — - — p)  cos  <p  u.] 
—2  pla  —  P  a  1  a  ^<p 


e.,  =  h„  x  u 

—1  -1  z 


_  ( 


(81) 


(82) 


(83) 


(84) 

(85) 


Its  modal  wave  number  is 


3.832,2,1/2 


rar.\ 


and  its  modal  characteristic  admittance  is 


a>£ 

Y  =  - - 

2  k„ 


(87) 


Other  parameters  of  interest  are  the  modal  voltages  V  ^  (0)  and 
Substituting  (80)  and  (84)  in  (23) ,  one  obtains 


V  ,  ,(0)  - 
111  a 


(88) 


^ 


(89) 


Now,  the  turns  ratios  n^  and  n^  are  readily  determined  from  (88)  and 
(89)  by  (78).  We  have 

0.816 

n  =  - 

11  a 


(90) 


_  0-991 
'l2  a 


(91) 
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VIII.  RESULTS  AND  CONCLUSIONS 


The  transfer  admittance  as  given  by  (34)  and  the  transmission 
cross  section  as  given  by  (63)  for  normal  incidence  and  one  propagating 
mode  in  the  cylindrical  cavity  were  calculated  and  selected  results 
are  presented  and  discussed  below.  In  the  following  explanations  of 
results,  an  effort  is  made  to  outline  how  quantities  of  interest  depend 
upon  electrical  dimensions  of  the  aperture-cavity-aperture  problem,  to 
examine  the  effects  of  different  dielectrics  in  the  cavity  region  b 
and  in  the  half-space  region  c,  and  to  look  at  the  effects  of  lossy 
material  filling  the  cavity.  Also,  the  influence  of  oblique  incidence 
(TE  field)  is  mentioned.  Finally,  the  case  of  two  propagating  modes 
in  the  cavity  is  considered  and  examined  by  a  numerical  example. 

The  case  of  normal  incidence  and  one  propagating  mode  in  the 
cylindrical  cavity  is  treated  in  Figs.  11-16.  The  parameters  Y^»  the 
transfer  admittance,  and  TA,  the  transmission  cross  section,  are  periodic 
functions  in  d,  the  depth  of  the  cavity.  Thus,  for  practical  consideration, 
they  arc  only  plotted  in  the  neighborhood  of  the  first  resonant  thick¬ 
ness.  The  permeability  of  all  regions  is  that  of  free  space  and  the 
permittivity  of  the  different  regions  is  specified  for  the  different 
cases.  The  apertures  A^  and  A^  are  assumed  to  be  identical  (i.e.  R^  =  R£) . 
Finally,  from  (40)  and  (86)  one  readily  obtains  that  the  TE^  mode  is  ‘the 
only  propagating  mode  provided 
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1.841 

2n 


1/2 

X 


< 


The  real  and  imaginary  parts  of  Y  ^ 
0.5  X^  and  for  aperture  radii  of  0.03 


a  < 


3.832 
2t r 


e  1/2 


X 

a 


are  shown  in  Fig.  11  for 

X  ,  0.05  X  ,  and  0. 1  X  , 
a  a  a 


where 


(92) 


0.03  X 
0.05  X 


various  radii  R,  and  R0.  Cases  shown  are 


A  is  Che  wavelength  in  the  half-space  region  a.  As  expected  from  (36) 

cl 

and  (37)  j  the  imaginary  part  of  becomes  dominant  as  the  apertures 
become  electrically  smaller. 

The  transmission  cross  sections  (63)  for  the  same  problem  are  shown  In 

2 

Fig.  12.  For  electrically  small  apertures  TA/X^  becomes  a  maximum,  as 

predicted  by  (61),  at  so-called  resonant  depths,  which  approach  A^/2,  where 

A,  =  1.234A  is  ttie  wavelength  of  the  TE. .  mode  for  £,  =  £  ,  V,  =  vj  ,  and 
La  J.-A  d  a  □  a 

2 

a  =  0.5A  .  The  peak  value  of  TA/X  for  all  three  cases  is  10  log  (3/4Tr)  = 
a  a 

-6.221  dB,  which  is  the  resonant  result  implied  by  (62). 

When  the  cavity  is  loaded  with  different  dielectrics  the  transmission 

resonances  occur  at  cavity  depths  farther  from  A^/2,  as  given  by  (61). 

This  fact  is  demonstrated  for  a  =  0.32A  ,  R.  =  R„  =  0.05X  and  £  -  c  =  £ 

a  1  2  a  a  c  o 

and  =  £  ,  2co, and  3e  in  Fig.  13a,  Fig.  13b  and  Fig.  13c,  respectively. 

Note  that  A  =  2.489)  0.928A  .and  0.680A  for  the  cases  of  e,  =  c  ,  2e  .and 

1  a  a  a  boo 

3t:  and  thus  the  resonance  actually  occurs  at  different  cavity  depths 

measured  in  wavelengths  in  the  half-space  region  a.  The  choice  of  the 

cylindrical  cavity  radius  a  =  0. 3° X  ensures  one  propagating  mode  in 

all  three  cases.  As  the  material  filling  the  cavity  becomes  more  dense,  the 

transmission  resonances  become  wider.  For  instance,  the  variation  in  d  for 

which  the  transmission  cross  section  is  not  loss  than  -60  dB  are 

0.996  d  <  d  <  1.004  d  ,  0.990  d  <  d  <  1.010  d  ,  and  0.987  d  <  d<  1.014 
res  res  res  res  res 

for  i ^  =  e  ,  t'k  =  2 and  =  3£o,  respectively.  Note  that  the  magnitudes 

of  the  peaks,  however,  are  constant.  This  result  is,  of  course,  in 

agreement  with  (58)  which  is  independent  of  C  . 

b 
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The  effects  of  a  different  dielectric  filling  the  half-space  region  c 


are  depicted  in  Fig.  14.  Here  the  cavity  radius  is  a  =  0.5Aa>and  the  radii 
of  the  apertures  are  R^  =  =  0.05Aa.  As  expected  from  (59), the  position 

of  the  resonance  condition  is  not  a  function  of  £  ,  and  therefore  is 

c 

not  affected.  The  peak  values  of  the  transmission  cross  section,  however, 
decay  in  magnitude  as  becomes  larger.  The  result  is  in  agreement 
with  (57)  and  (58)  which  reduce  here  (y^  =*  y  =  1  and  x  =  1)  to 


Re(Y  )  =  (-1)®  +  1  +  E  3/2) 

12  res  ?  c 

3n  A^ 
a  a 


T 

res 


3/2  X2 

[1  +C£  3/2lr^ 

c 


(94) 


where  £  =  £  /e  as  defined  by  (53). 

c  c  a 

If  the  dielectric  filling  the  cavity  is  lossy, the  expected  decrease 

in  the  transmission  cross  section  peaks  is  seen  in  Fig.  15,  which  presents 

the  lossless  case  £,  =  £  ,  the  lossy  case  £.  =  (1  -  i0.001)£  ,  that  is, 
bo  o  o 

with  Q  factor  =  1000, and  the  lossy  case  £^  =  (1  -  j0.01)£o>  that  is,  with 

Q  factor  =  100.  Here  the  cavity  radius  is  a  =  0.5Aa.  Note  that  the  loss 

does  not  affect  the  position  of  the  resonances.  The  real  and  imaginary 

parts  of  for  this  case  are  shown  in  Fig.  16.  As  may  be  observed,  the 

loss  does  not  affect  the  imaginary  part  of  Y.^*  In  contrast,  its  real 

2 

part  has  changed  significantly.  Re(Afl  Y12>  for  the  case  £^  =  (1  -  j0.01)£o 
could  not  be  shown  in  Fig.  16a  since  it  varies  between  19  and  29  mhos. 

If  the  two  apertures  are  of  different  size,i.e.,  4  R2,  (57)  and 
(58),  for  the  case  of  three  similar  media  (i.e.,  =  yc  =  £c  =  1), 

reduce,  respectively,  to 

Rett,,)  -  (-Dm  (2  -  (1  +  x)d  +  £)]  (95) 

12  res  3n  yf-  X 

a  a 
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Fig.  15.  Cases  shown  are 


7T[2  -  (1  +  x)(l  +  ^)]2  A 

From  (96)  one  can  easily  show  that  Trgs  is  maximum  when  y  =  1,  that  is 
when  the  two  apertures  are  identical.  This  is  in  agreement  with  the  fact 
that  when  R^  =  the  power  transmitted,  at  resonance,  is  equal  to  the 
maximum  power  that  can  be  delivered  to  a  load  by  the  source  and  its  internal 
admittance  +  jB  ^  in  Fig.  7. 

The  effect  of  TF.  (transverse  electric  to  the  cylinder  axis)  oblique 

incidence  is  shifting  the  transmission  cross  section  graphs  down  by 
2 

10  log(cos  0jnfi)  decibels,  where  8^  ,  the  angle  of  incidence,  is  shown  in 

Fig.  8.  Otherwise,  TF.  oblique  incidence  is  the  same  as  TEM  normal  incidence. 
Hence,  no  further  detail  regarding  TE  oblique  incidence  is  necessary. 

Finally,  the  case  of  normal  incidence  and  two  propagating  modes  in 

the  cylindrical  cavity  is  treated  in  Fig.  17.  The  permeability  and 

permittivity  of  all  regions  are  those  of  free  space.  The  apertures 

and  A  are  identical  and  R  =  R„  =  0.05X  .  From  (86)  and  the  fact 
L  12  a 

that  the  wave  number  of  the  third  mode,  the  TE^  mode,  is 


,  r  2  ,5. 331, 2,1/2 

k3  =  [<,  ubeb  -  ( — t; — )  I 

one  readily  obtains  that  the  TE^  and  TM^  modes  are  the  only  propagating 
modes  provided 

3  832  Ga  5  331  ^a 

tP  <i?>  V  a  <  x, 

b  b 


Accordingly  ,  the  cavity  radius  is  taken  to  be  a  =  0.7SX  . 

a 

The  transmission  cross  section  in  Fig.  17  is  plotted  in  the  neighbor¬ 
hood  of  the  first  two  resonant  depths.  The  transmission  cross  section 
peaks  are  located  at  these  depths,  which  approach,  respectively,  X  /2 
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Fig.  17.  Plots  of  transmission  cross  section  versus  d/X  for  e  =  e  =  c  = 

a  a  b  c 

a  =  0.75Xa>  and  =  0.05Xa  for  the  case  of  two  propagating  modes 

in  the  cavity. 


and  \^/2,  whore  A^  =  l.086A^  is  the  wavelength  of  the  TE^  mode  and 

A  =  1.718A  is  the  wavelength  of  the  TM  mode.  The  magnitude  of  the 
Z  3  11 

two  peaks  is  the  same  and  equal  to  -6.221  dB,  which  is  the  peak  value 
obtained  when  only  one  mode  was  propagating. 


IX.  DISCUSSION 


Electromagnetic  coupling  between  two  half-space  regions  through 
an  aperture  to  cavity  to  aperture  system  has  been  investigated  and  a 
simple  equivalent  circuit  has  been  developed  for  the  special  case  of 
a  circular  cylindrical  cavity  with  a  small  circular  aperture  centered 
in  each  of  its  end  faces.  It  was  found  that  for  certain  cavity  depths 
the  system  becomes  resonant,  thereby  increasing  by  orders  of  magnitude 
the  power  transmitted  by  the  aperture  to  cavity  to  aperture  system  over 
what  it  would  transmit  for  other  cavity  depths. 

It  should  be  emphasized  that,  if  the  generalized  concept  of  small- 
aperture  polarizability  which  adds  a  radiation  conductance  term  to  the 
aperture  admittance  had  not  been  used,  infinities  in  the  aperture  fields 
and  power  transmitted  would  have  occurred . 

The  result  that,  at  resonances , the  transmission  cross  section  of  the 
coupling  system  is  independent  of  the  size  of  the  apertures  (provided 
R1  =  R2^  *S  to  ke  exPecte<i  because  similar  results  have  been  obtained 
for  other  transmission  and  scattering  problems.  For  example,  the  trans¬ 
mission  width  of  a  narrow  infinitely  long  slot  in  a  thick  conductor  at 
resonance  is  ^  ,  regardless  of  its  actual  width  [A],  Moreover,  the 
transmission  cross  section  obtained  here  is  the  same  as  that  obtained 
in  [7]  for  a  small  slot  in  a  zero-thickness  plane  conducting  screen 
resonated  by  a  capacitor  placed  across  its  midpoint. 

The  preceding  paragraph  applies  only  to  ideal  loss-free  problems. 
Dielectric  and  conductor  losses  can  significantly  decrease  the  trans¬ 
mission  cross  section  at  resonance.  This  decrease  has  been  demonstrated 
via  numerical  examples. 
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